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Abstract The maximally extended Reissner-Nordstrom (RN) manifold with e 2 < m 2 begs for attaching a 
material source to it that would preserve the infinite chain of asymptotically flat regions and evolve through 
the wormhole between the RN singularities. So far, the attempts were discouraging. Here we try one more 
possible source - a solution found by Ruban in 1972 that is a charged generalisation of an inhomogeneous 
Kantowski-Sachs-type dust solution, ft can be matched to the RN solution, and the matching surface must 
stay all the time between the two RN event horizons. However, shell crossings do not allow even half a cycle 
of oscillation between the maximal and the minimal size. 



1 Motivation 



Avoiding the Big Bang singularity in cosmological models had been a recurring idea in the literature. The 
hopes for constructing a model without singularity were largely dashed by the singularity theorems of 



Hawking and Penrose (see Ref. [T] for a review). They were temporarily revived by the finding of Vickers 
[2] that in a charged dust model generalising that of Lemaitre [3] and Tolman [3] (LT) the Big Bang 
can be prevented by the charge distribution^] provided that the absolute value of the charge density is, in 



> 

qv ' geometrical units, smaller than the mass density. This finding created another expectation, that a finite ball 

of charged dust matched to the Reissner [U] - Nordstrom [7 (RN) solution would be able to collapse and 
bounce through the wormhole of the maximally extended RN manifold, thereby giving physical meaning to 
the infinite chain of black holes and asymptotically flat regions) 2 ] However, this renewed hope was dashed 
again in 199f by Ori [12] who proved that, under exactly the same conditions that prevent the Big Bang, shell 
crossings will inevitably appear and destroy the dust ball before it enters the wormhole. Then, Krasinski 
and Bolejko [13,14 found a gap in Ori's assumptions and tried to improve upon his result. Namely, Ori 
assumed that the ratio of the absolute value of the charge density to the mass density (call this ratio a) 
is smaller than 1 everywhere, including the centre of symmetry. Krasinski and Bolejko considered the case 

p^ ' when a — > 1 at the centre, while being < 1 everywhere else. The situation turned out to be better, but not 

definitively. With initial conditions carefully tuned, the charged dust ball could go through the wormhole 
just once, being destroyed by shell crossings soon after reaching the maximal size in the next asymptotically 
flat region. Moreover, a direction-dependent singularity necessarily appeared at the centre of the ball at 
the instant of minimal size. This is not a satisfactory situation, but the best result achieved so far. 
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1 This charged generalisation of the LT model was first found as a solution of the Einstein - Maxwell equations by 
Markov and Frolov in 1970 [5]- 

2 The extension of the RN solution composed of two coordinate patches was first calculated by Graves and Brill in 1960 
PJ]. The infinite mosaic of conformal diagrams shown in Fig. ^ first appeared in the paper by Carter |9] in 1966, and was 
described in detail in a review article by Carter 1101 . See Ref. Ill| for another pedagogical introduction. 



In an attempt to find a better model, we now tried to match the Ruban [15] charged dust solution 
to the RN metric and see what results. The charged Ruban solution is a generalisation to nonzero charge 
of the nonstatic dust solution investigated earlier by Ruban [16U17] , which in turn is an inhomogeneous 
generalisation of the Kantowski - Sachs model [TH]o The matching of these two solutions is very easily 
achieved, with an interesting result. The Ruban solution represents a pulsating charged dust ball whose 
outer surface must remain forever between the two event horizons of the RN solution. It touches the inner 
event horizon at its minimal size and the outer horizon at its maximal size. However, a simple investigation 
shows that shell crossings are inevitable inside the ball and do not allow even a half cycle of an oscillation, 
appearing between the maximum and minimum size both in the expansion phase and in the collapse phase. 

Thus, the field is still open for finding a physical nonstatic source for the RN solution that would 
proceed through the wormhole. The next thing to do is to find a charged perfect fluid solution of the 
Einstein - Maxwell equations in which pressure gradients would prevent the formation of shell crossings. 
This, however, promises to be an extremely difficult task. 

In this paper, we first present the RN solution in coordinates adapted to the matching, then the charged 
Ruban solution as found by the author, then we prove that the matching can be done, and finally we prove 
that shell crossings in the Ruban solution are inevitable. We also investigate the limit of zero charge of 
the Ruban - RN configuration. Then the dust source is matched to the Schwarzschild solution across a 
hypersurface that remains inside the Schwarzschild horizon, only touching it from inside at the moment of 
maximal expansion. Shell crossings can then be avoided by an appropriate choice of the arbitrary functions 
in the Ruban model. However, this model has a finite time of existence, being born in a Big Bang that 
matches to the past Kruskal [21 - Szekeres [22] singularity and crushing into a Big Crunch that matches 
to the future KS singularity. 

We use the signature (H ). We do not assume "units in which G = c = 1" ; whenever these constants 

are absent, this means they were absorbed into other symbols. For example, our time coordinate will be 
t = ex [the physical time], our M will denote (G/c 2 )x [the mass], etc. The labelling of the coordinates is 
{x°,x 1 ,x 2 ,x' i ) = (t,r,$,<p). 

2 The Reissner— Nordstrom solution with e < m and its maximal extension 

The Reissner [6] - Nordstrom [7] solution in its standard form is 

ds 2 = Fdt 2 - (1/F)dr 2 - r 2 (dtf 2 + sin 2 -ddtp 2 ) (2.1) 

where 

„, . def _. 2m e 2 1 2 , n „> 

F(r) = 1 + -2 + 0^ 2.2 

r r z 6 

(we include A for a while). For reference, Fig. Q] (adapted from Ref. |11|) shows the Carter - Penrose diagram 

of the maximal analytic extension of the underlying spacetime for the case A = and e < m , to which we 

will limit our attention in the main part of this paper. The ragged lines in the diagram are the curvature 

singularities at r = 0, the lines marked "inf" are the past and future null infinities, r+ are the outer event 

horizons at 

r = r+ = 771 + v m 2 — e 2 , (2-3) 

and r_ are the inner event horizons at 



r = r_ = m — \/m 2 — e 2 . (2-4) 

These two values of r are zeros of the function F when A = 0. Regions I and III are asymptotically flat, 
regions II and IV are contained between the two horizons. The thin curved lines are those on which r 
is constant; they are timelike in the asymptotically flat regions where r > r+ and inside the inner event 
horizon where r < r_; between the horizons where r_ < r < r+ they are spacelike. 

In the following we will be interested in the region between the horizons, where F < 0, so r in (|2.1[) 
becomes the time coordinate and t becomes a spacelike coordinate. It is easy to verify that the lines on 
which (t,-ff,ip) are all constant are geodesies. 

For later convenience, we rename the coordinates in (12.111 as follows 



(t,r) = (p,R), (2.5) 



3 The neutral dust solution investigated by Ruban was first found as a solution of the Einstein equations by Datt in 1938 
|19l . but instantly dismissed as being of "little physical significance". 




Fig. 1 The conformal diagram of the maximally extended Reissner-Nordstrom spacetime with A = and e 2 < m? . It 
can either be interpreted as an infinite chain of copies of the segment within the rectangle, or one can identify the upper 
side of the rectangle with the lower side, obtaining a manifold with closed timelike and null lines. More explanation in the 
text. 



and then transform the il-coordinate in the region F < to r(R) defined by 

dR Y - F - ^ + 2m e " l AR 2 
d7j — F — 1 + ^-W-3 AR - 

After this, the metric in this region becomes 

ds 2 = dr 2 - [-F(R(r))]dp 2 - R 2 (r) (dtf 2 + sin 2 tidtp 2 



with R(t) being defined by (f2~ 

When A = 0, eq. (|2.6[) can be integrated to give the following explicit function t(R): 



t-tq = n 



1 



di? 



2m e 2 

1 + 1T-W 



i\/—R 2 + ImR — e 2 + fmi arcsin 



R — rn 



\f~ti 



(2.6) 



(2.7) 



(2.8) 



where fi = ±1 and tq is an arbitrary constant. 

For some purposes it is more convenient to introduce the parameter r/ by (R — rn)/\/m 2 — e 2 = sin(?i — 
7r/2) = — cosr?, and then 

R = 771 — \/m 2 — e 2 cos i), t — to = /x I ra?i — v m 2 — e 2 sin 77 1 , (2-9) 

where to = to — fimn/2. This shows that as r increases, _R is oscillating between the r_ and r+ given by 
dO and J2" 



3 The charged Ruban solution [15J 

To derive this solution, we assume comoving coordinates, spherical symmetry and a charged dust source. 
For the metric we assume a less general form than spherical symmetry would allow, namely 

ds 2 = e C(t ' r) dt 2 - e A(t ^dr 2 - R 2 (t) [<M 2 + sin»d^ 2 ] , (3.1) 

where C(t,r), A(t,r) and R(t) are functions to be found from the Einstein - Maxwell equations. The 
limitation of generality is in R(t); in the most general case it would depend on r as well and the field 
equations would lead to a charged generalisation of the LT model. 

Details of the calculation are given in Ref. [11 . The field equations show that R,t C, r = 0. With R,t = 0, 
we obtain an (electro-) vacuum solution, which is presented in Appendix [X]- it is a coordinate transform 
of the well-known Robinson solution [20]. Thus C, r = 0, which means that the dust is moving on geodesies. 
A transformation of t can then be used to achieve C = 0. The other equations imply 

*rf a = - 1+ ir-i* ~I ar2 > (3 - 2) 

where M is a constant of integration, and Q is the electric charge within the sphere of coordinate radius r. 
Since all other quantities in this equation are independent of r, it follows that Q must be constant. This 
implies that the charge density is zero everywhere. Thus, the dust particles must be neutral, they only 
move in an exterior electric field. 

For the other metric function we obtain the following solution: 



e A/2 = R„ 



M "/557 + y(r) 



(3.3) 



where X(r) and Y(r) are arbitrary functions, and the expression for matter density in energy units is 

2X 

R 2 e A/2 



(3.4) 



Equations (|3.1[) - (|3.4p define the charged Ruban solution, first semi-published in 1972 [15], and then 
mentioned in a later paper [23] ■ Note that (|3.2[) is identical to (|2.6[l . so the function R(t) defined by (|3.2[) 
will be the same as R(t) defined by (|2.6(1 . Note also the following: 

1. With X = 0, eq. (I3.4|) shows that the Ruban solution becomes electro- vacuum; in fact it then becomes 
the Reissner - Nordstrom solution expressed in the (r, p) coordinates defined by (|2.5|) - (|2.6[) . 

2. When Y = BX, where B is a constant (B = being allowed), the r-dependence in (|3.4|) cancels out 
and by defining a new r by r = J X(r)dr we make also the metric independent of r. The spacetime then 
becomes spatially homogeneous with the Kantowski - Sachs symmetry; in fact it is then the generalisation 
of the Kantowski - Sachs solution to nonzero charge and cosmological constant. 

3. The geometry of a 3-space t = const in (|3.1[l is that of a 3-dimensional cylinder whose sections r = 
const are spheres, all of the same radius, and the coordinate r measures the position along the generator. 
The space is inhomogeneous along the r-direction, and the electric field has its only component also in the 
r-direction. The radius R of the cylinder evolves with time according to Eq. (|3.2p . 

4. The matter density in this solution, given by (|3.4p . depends on r and is everywhere positive if X > 0. 
Thus, the amount of rest mass contained inside a sphere r = ro = const does depend on the value of ro, and 
is an increasing function of r. Nevertheless, as seen from f|3.2[> . the active gravitational mass M that drives 
the evolution is constant. Ruban [T7] interpreted this property as follows: the gravitational mass defect of 
any matter added exactly cancels its contribution to the active mass. 



4 Matching the Ruban and RN solutions 



To match two spacetime regions to each other we have to prove that on the hypersurface H forming the 
border between them both 4-dimensional metrics induce the same 3-dimensional metric and the same second 
fundamental form. The formulae we use here are derived in Ref. |11| . 

We will show that the RN metric in the coordinates of (|2.6p - (|2.7p and the Ruban metric defined by 
(|3.1[l - (|3.3p with C = can be matched along any H of constant p = r. On each such H the 3-metrics will 
be identical when r and i are identified and 



m = M, e = Q, 

since then the R(t) of (|2.7p and the R(t) of (|3.ip will obey the same equation, 
The coincidence of the second fundamental forms requires that at H [11] 



1 dhi 



N dr 



1 dhi 



N dp 



(4.1) 
2.6p and (|3.2p respectively. 



(4.2) 



where i,j = 0,2,3 and M is the r-component of the unit normal vector to H , thus \g r r\ A/" = 1 on each side 
of H. But the relevant components of the Ruban metric do not depend on r, and the relevant components 
of the RN metric do not depend on p, so (|4.2|) is fulfilled in a trivial way: it reduces to = 0. Note, however, 
that the matching is possible only in that region of the RN manifold where F < in (|2.ip - (|2.2p . In the 
subcase 4 = 0, this is the region between the two event horizons. 

Equation (|2.9p applies also in the Ruban region, with m = M and e = Q. It is independent of r, so 
each constant-r shell evolves by the same law. This means that all shells, including the outer surface of 
the Ruban region, oscillate between R = r_ and R = r + , where r_ and r+ are given by (|2.3p - (|2.4p . If 
the solution could be continued to infinite values of t, the conformal diagram of the complete manifold 
would look as in Fig. [2] However, shell crossings make the completion of even half a cycle of oscillations 
impossible, see next section. 



5 Shell crossings in the charged Ruban solution are inevitable 

For the rest of this article we assume 4 = 0. In this case the Ruban metric is 

ds 2 = dt 2 - e A(t ' r) dr 2 - R 2 {t) (dtf 2 + sin 2 dd V 2 ) , 

where the solution for R(t) can be written as (cf. eq. (|2.9p ): 

R = M - y/M 2 -Q 2 cosr], t - t = fi (Mr] - a/A/ 2 - Q 2 sin 77) , 



(5.1) 



(5.2) 



where p = +1 in the expansion phase and p = — 1 in the collapse phase. With 4 = 0, the function e ' 2 can 
be calculated explicitly. For this purpose, we find from (|3.2p 



(5.3) 



R,t=(i\ -1 + 



2M Q 2 
~R~^R^' 



Note that this solution exists only when Q 2 < M 2 O Now we use this in the integral in (|3.3p rewritten as 

follows 

dt f dR (54) 



RR.-, 



RR.-, 



Then we obtain 



^A/2 



X{r) 



n -2,1- MR I 2M Q 2 . R-M 

2 + Q 772 7v? - \ /- 1 + -FT - 1*2 arcsm / 

M 2 -Q 2 V R R 2 \ ,/M 2 - Q 2 



1 2M Q 2 



(5.5) 



4 With Q 2 = M 2 , Eq. 115.31 1 has the static solution R = M = \Q\. As remarked above A3. 211 . this leads to the Robinson 
solution 1201 ; see Appendix \K\ 




Fig. 2 The maximally extended Reissner-Nordstrom spacetime with a hypothetical charged matter source matched to 
it. The extent of the source in this diagram is unknown, even if it is the charged Ruban solution - because of the arbitrary 
functions it contains. This is why the left edge is schematically marked with a dashed line. The two long wavy curves show 
possible matchings at two different values of the radial coordinate; this coordinate is the r of 115. f II and the p of 112.711 . Note: 
these are the curves on which r and p are constant. It is R that pulsates between the r_ and r-f. of 112. 311 and 112.411 . and 
it changes along these curves. The lines of constant R are the horizontal hyperbola-like arches. If the matter source is the 
charged Ruban solution, then shell crossings appear within each half-cycle of oscillation. The positions of two of them are 
schematically marked with thin horizontal straight lines. 



A shell crossing is where e ' 2 = 0, because at such a location two matter shells whose r-coordinates 
differ by dr are at a zero distance, as seen from (|5.ip . and stick together. From (|3.4p one sees that at such 
a location the mass density becomes infinite, so this is a curvature singularity. The question we seek to 
answer is now: can the functions X{r) and Y(r) in (|5.5p be chosen in such a way that e ' 2 ^ everywhere. 

Unfortunately, the answer given by the formulae above is a definitive "no" , i.e. shell crossings are 
inevitable. To see this, let us first recall that R changes between the values 



R- = M - v^ 



and R+ = M+ ^ M 2 - Q 2 



(5.6) 



(this is a copy of (|2.3f) - (|2.4p rewritten in the notation for the Ruban model). At both these values the 
expression under the square root in (|5.5|) is zero, and in the whole range < R- < R < i?+ the function 
e ' 2 given by (|5.5() is continuous. Now, at R = R- we have 



e A/2 \ =-X(r)— F JL , (5.7) 



\R-R- 



VM 2 



while at R = R+ we have 



e A/2\ =X{r) R+ = (5g) 

\R=R+ ^/KP - Q 2 

Thus, whatever the sign of X(r), the signs of e ' 2 at the ends of the range [_R_,_R+] are opposite, i.e. 
e ' 2 = for some value of R within this range. This is a shell crossing. Note that it appears every time 
when R traverses this range, which means that the shell crossings will not allow R to go through even half 
a cycle of oscillation. 

The geometrical nature of this shell crossing is different than in the LT and Szekeres models [241125] , 
where they have been investigated quite thoroughly [26,27,28.. In the LT and Szekeres models, the spheres 
that collide are one within the other. At a shell crossing, the smaller sphere catches up with the larger one 
during expansion (or vice versa during collapse). In the Ruban model, the spheres are surfaces of constant r 
in the 3-dimensional cylinder t = constant, and they move up and down along the generators as the cylinder 
expands or collapses. It turns out they will collide before the cylinder manages to proceed all the way from 
the minimal radius to maximal, or from maximal to minimal. One more property of such a shell crossing 
is noteworthy: its locus does not depend on r, which means that all the spheres with different values of r 
collide at the same moment. 

Kantowski and Sachs [18] noted the possible occurrence of this singularity in the subcase Q = Y(r) = 
of the Ruban model, where they said "The singularities for the closed models are of two kinds. (...) In one 
kind of a singularity the cylinder squashes to a disk, in the other it contracts to a line." What we called 
shell crossing here is the "squashing to a disk" . 

The shell crossings could possibly be prevented by pressure with a nonzero gradient in the r-direction, 
but such exact solutions are not yet known. Then, assuming that the pressure and its gradient would be 
zero at the surface of the charged fluid ball, its surface would follow a timelike geodesic in the RN spacetime, 
and the conformal diagram could really look like Fig. [2] 



6 Absence of shell crossings in the Ruban solution with zero charge 

When Q = 0, the charged Ruban solution goes over into the Datt - Ruban neutral dust solution [19U171 
[TTj . Then from (|2.9p , adapted to the notation for the Ruban solution, we obtain 

R = M(l-cosr)), t-to = fxM(ji-smrj). (6.1) 

This shows that now R starts from at t = to, then increases to 2M at t = to + fxMw, and decreases to 
again at t = to + 2(iMtt. At R = the model has a Big Bang/Crunch type singularity. Equation (|5.5|) 
simplifies tqf| 



e^ = ,/-l + ^ JTM 



.1 -1- £M- 

1 T R 



arcsm -—- — 1 
2M V M 



+ „Vi,>}> ^ \- 1+ ™jr(fi,r). (6.2) 



Now the spacetime is singular at R = 0, so the model exists only for the finite time T = 2-rvM between 
the Big Bang and Big Crunch, but, unlike in the charged case, the functions X(r) and Y(r) can be chosen 
so that within this interval there are no shell crossings. We show below that, with an appropriate choice, 
T > during the whole evolution. 

Namely, with R = at the start of the expansion phase, where /x = +1, we have 

T\ R = = X^ + Y (6.3) 



5 Equation (16.2ft is equivalent to (19.101) in Ref. |11| under the renaming Y = fj,(Y — 7r/2), in virtue of the identity 
arcsin (R/M - 1) + tt/2 = 2 arcsin y/R/(2M). 



(6.4) 



Thus, X > and Y > will guarantee that F\ r= q u =+i -* 0- We also have 

dJ" _ X 

which will be positive for all < R < 2M if X > 0, becoming +oo as R — > 2M. Since also J" — ► +oo as 
i? — > 2M, the switch from expansion to collapse, i.e. from /i = +1 to /i = —1, will keep J" positive at the 
beginning of collapse. So we only have to guarantee that T > at the end of the collapse phase, when 
R — > again. This will be the case when 

T\ R=0 ^ = _ X = X^-Y>0. (6.5) 

Consequently, both (|6.3p and (|6.5p will be positive during the whole cycle when 

0<F<X| (6.6) 

at all values of r. The fact that the shell crossings become absent when Q = shows that the limiting 
transition Q — > is discontinuous; just as it was in the RN — >• Schwarzschild transition. 

The neutral Ruban solution can be matched to the Schwarzschild solution, as shown by Ruban himself 
[17] . The matching hypersurface stays within the Schwarzschild radius except at R = 2M, when it touches 
the horizon from inside. However, this configuration has a finite time of existence, just as the Kruskal [21] 
- Szekeres [22 manifold between its two singularities. 

7 Summary 

We investigated the charged Ruban solution [15] as a possible matter source for the maximally extended 
Reissner - Nordstrom solution with e 2 < m 2 . The matching of these two solutions is easily achieved. The 
hypersurface that forms the boundary between the Ruban and RN regions stays all the time between the 
two RN event horizons, touching the inner one at its minimal radius and the outer one at maximal radius. 
However, shell crossings make the completion of even half a cycle of such an oscillation impossible; they 
appear between each [r_,r-|-] and each [r_|_,r_] pair. 

Shell crossings can be avoided in the limit of zero charge, when the arbitrary functions in the Ruban 
solution obey a simple inequality. Then the Ruban solution can be matched to the Schwarzschild solution 
inside the event horizon, but the model has a finite time of existence. 

The problem of finding a matter source to the maximally extended RN solution thus still remains open, 
with one more possibility being now elliminated. 

A The limiting case Q 2 = M 2 , A = 

As stated in the footnote below 115.31 1. this leads to the static solution R = M = \Q\. This limit is not a subcase of 13.311 . 
and the Einstein - Maxwell equations have to be solved anew. With R,t = they lead to the electro-vacuum solution, in 
which 

e A / 2 = a(r) cosftt/Af) + 0(r) sin(t/M), (A.l) 

where a(r) and j3(r) are arbitrary functions, and the metric is 

ds 2 = dt 2 - e A dr 2 - A/ 2 (dtf 2 + sin 2 i?d^ 2 ) (A. 2) 

(we recall that M = \Q\, where Q is the electric charge). This should be compared to the Nariai solution |29 | in the form 
found by Krasifiski and Plebanski |30| : 

ds 2 = [a(t) cos(p/£) + b(t) sin(p/£)] 2 di 2 - dp 2 - £ 2 (dtf 2 + sin 2 tfd>p 2 ) , (A.3) 

where £ is an arbitrary nonzero constant. This is a vacuum solution with the cosmological constant A = l/£ 2 . Equation 
1A.3I I is related to llA.lll - llA.211 by the interchange of t and r, and has the same geometrical structure - it is a Cartesian 
product of two surfaces of constant curvature. 

In fact, HA. Ill — llA.21 1 is a coordinate transform of the Robinson solution 20 . Namely, by a rather complicated sequence 
of coordinate transformations, strictly analogous to those used in Ref. 1301 for the Nariai solution, I IA.U — (IA.2I ) may be 
transformed to 

£ 2 
ds 2 = —2 \dl- 



dr' 2 - r' 2 (dtf 2 + sin 2 tfdip 2 )] , (A.4) 

which is identical with one of Robinson's formulae, with Robinson's A = l/£. 
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